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Models of the iron Ka fluorescence line and the Compton 
shoulder in irradiated accretion disk spectra 
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ABSTRACT 

We present a full set of model atmosphere equations for the accretion disk around a 
supermassive black hole irradiated by hard X-ray lamp of power-law spectral distribu- 
tion. Model equations allow for multiple Compton scattering of radiation on free elec- 
trons, and for large relative photon-electron energy exchange at the time of scattering. 
We present spectra in specific intensities integrated over the disk surface. Theoretical 
outgoing intensity spectra show soft X-ray excess below 1 keV, and distinct Ka and 
Kf^ fluorescent lines of iron. We demonstrate the existence of the Compton Shoulder 
and claim that it can contribute to the asymmetry and equivalent widths of some 
observed Fe Ka lines in AGN. Our models exhibit the effect of limb-brightening in 
reflected X-rays. 
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■ 1 INTRODUCTION 

Observations of many active galactic nuclei (AGN) directly 
I show the existence of an X-ray source in the innermost re- 
. gion of AGN, which emits radiation with power-law spec- 
' trum extending to above 100 keV (|Guainazzi et al.l 1 19991 : 
, iDone et al.l |2003| ) . The exact location of the source is still 
' unclear, nevertheless, such hard radiation strongly affects 
' other emitting regions including nearby accretion disk. Ef- 
fects of the redistribution of hard X-rays into the soft X-ray 
and ultraviolet domains are frequently observed in many 
AGN startin g with the paper Pounds et al.l (|l990t ) up to the 
recent paper iNandra fc Iwasawa l|2007l ). 

The most important spectral feature indicating that the 
reflection occurs from an accretion disk, is a fluorescent iron 
Ka line. Close to the supermassive black hole, the emis- 
sion iron line profile is relativistically broadened and skewed 
l|Fabian et al.lll989l : iRevnolds fc Begelmanlll997l ). which to- 
gether with a circular gas movement in a disk, makes final 
profile asymmetric with a strong red wing. 

The detection of iron line profile was a goal of sev- 
eral X-ray satellites, and the results depend on the res- 
olution of a particular instrument. The best known ob- 
ject with broad iron Kg line, MCG-6-30-15 was o bserved 
by ASCA ("Tanak a et all 1 19951 : llwasawa et"all ll999h. In the 
Chandra and XMM era, line became narrow or contained 
both broad and narrow component (jReeves et al ] 12004 
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IPetrucci et"alll2007l ). Recently, more often authors are re- 
porting a detection of "iron line complex" additionally af- 
fected by intrinsic warm absorber fYagoob fc PadmanabhanI 
^2004; Risaliti ct al. 2005; Matt ot al. 2006). 

The most recent data coming from Suzaku look 
very promising IIMiniutti et al ] l2007l : iReeves et all l2007l : 
lYaaoob et aP l2007l l. The iron line complex contains at 
least two components: narrow Gaussian at 6.4 keV, and 
broad red wing, the last one usually fitted with a disk 
li ke shape, but also with Compton shoulder in the work 
of lReeves et al.l (|2007l ). Additionally iro n Kg line at energy 
7.057 keV is clearly seen in Suzaku AGN |Reeves et al ] |2007l : 
iMarkowitz et al.ll2007l : [Yaqoob et al.ll2007l b 

In this paper we present quite consistent and numer- 
ically exact modeling of the iron line complex originating 
from an accretion disk atmosphere irradiated by hard X- 
rays with power-law spectral distribution. Our computations 
of the vertical structure and outgoing spectrum of the disk 
are based o n rigorous m ethod of the theory of stellar at- 
mospheres l|Madeil [l99il : [Madei fc Rozariskal boOOal lbl). All 
existing opacity sources both true absorption and Compton 
scattering ar e treated equivalent l y in th e single code fully 
described in iMade i fc Rozaiiska' (|2004| . hereafter MR04). 
Adopting atmospheric computations for an accretion disk 
around supermassive black hole, we solve the structure and 
obtain outgoing spectra for eight neighboring rings. The fi- 
nal spectrum presented in specific intensity scale for differ- 
ent aspect angles is integrated over radii and presented at 
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the source frame. No kinematic special relativity effects are 
included in our model. 

Our radiative transfer equation includes effects of multi- 
ple Compton scattering of radiation on free electrons in rela- 
tivistic thermal motion, and rich set of bound-free and free- 
free opacities. We allow for a large relative photon-electron 
energy exchange at the time of Compton scattering and, 
therefore, are able to reconstruct Compton scattering of pho- 
tons with energy approaching or even exceeding the electron 
rest mass. In this paper the opacities are supplemented by 
formulae describing emission of fluorescent Ka and lines 
of low ionized iron. 

The structure of the paper is as follows: Sec.[2]describes 
sample accretion disk model, and Sec. [3] presents most im- 
portant equations used in our code. Results are presented in 
Sec. H] and Sec. [5] and summarized in Sec. |6] and Sec. [T] 



2 MODEL OF A SAMPLE ACCRETION DISK 

We consider an exemplary model of the accretion disk 
around supermassive black hole located at the center of 
AGN. Mass of the black hole equals Mbh = lO^'M©. We 
divide our disk onto 8 concentric rings situated at differ- 
ent distances from the center. We assume constant moder- 
ate value of accretion rate in each ring, equal rh — 0.03 in 
units of Eddington accretion rate with accreting efficiency 
7] = 1/12 suitable for the Schwarzschild black hole. 

At the first step, at each radius, we computed vertical 
str ucture of non-ir r adiate d disk using the method described 
in iRozanska et al.l (|l999l ). We integrated equations of the 
disk vertical structure (equation of the hydrostatic equilib- 
rium, equation of state, equation of local energy generation, 
transfer in diffusion approximation) assuming that viscosity 
is proportional to the total pressure, i.e. Ptot = Prad + Pgas- 
As the result, for each ring, we derived effective tempera- 
ture and vertical gravity which affects atmosphere, and we 
adopted them for actual computations. All disk parameters 
used in further computations are summarized in Tab. [1] col- 
umn 3 and 4. 

At the next step, we assume that the disk is irradiated 
by an external point like X-ray lamp, located above the first 
innermost ring, ri — 3.48rschw, at the height hi = drschw 
Irradiation spectrum is in the form of power-law with expo- 
nential low and high energy cut-off (see Sec. I3.1.2|) . 

The code bv lRozariska et al] l|l999l ) predicted the value 
of gravity, loggm, for each ring of the non-irradiated disk. 
When external irradiation was applied, the atmosphere al- 
ways lost its hydrostatic equilibrium. We had to increase the 
value of surface gravity by a factor of two or more to restore 
hydrostatic equilibrium in those rings, where equilibrium 
probably may not exist. There is no firm physical justifi- 
cation for such a step, of course. Values of gravity adopted 
for the advanced radiative transfer computations, log girr, 
are listed in the fifth column of Tab. [T] 

We take into account three different chemical compo- 
sitions of the accretion disk presented in Table [2] The first 
composition, Comp. I, consists only of hydrogen, helium and 
iron in proportions close to solar values. Comp. II is just the 
solar chemical composition for only ten most abundant ele- 
ments given bv IChaisson fc McMiUai] l|l99d) . We also con- 



Table 1. Parameters of the disk for Mbh = IO^'Mq, m = 0.03. 



ring No. ^'/i^Schw '^cff 

log gm log girr 



1 3.478E-I-00 1.018E-I-05 5.230E-I-00 5.5300 

2 4.134E-(-00 1.073E-f05 5.274E-I-00 5.5740 

3 5.204E-(-00 1.020E-f05 5.170E-I-00 5.4700 

4 6.552E-I-00 9.238E-I-04 5.001E-I-00 5.3010 

5 8.248E-I-00 8.181E-t-04 4.801E-I-00 5.1010 

6 1.038E-I-01 7.152E-f04 4.582E-(-00 5.1820 

7 1.385E-I-01 5.976E-f04 4.294E-(-00 5.594 

8 1.847E-I-01 4.950E-I-04 3.996E-I-00 5.596 



Table 2. Three chemical compositions used in our calculations. 
Entries of the table display number abundances of elements rela- 
tive to hydrogen. Three rows at the bottom present mass abun- 
dances X, Y, and Z for all three composition. 





Comp. I 


Comp. II 


Comp. Ill 


H 


1.0 


1.0 


1.0 


He 


1.1 X 10-^ 


9.54 X 10-^ 


9.54 X 10-2 


C 




4.72 X 10"* 


4.72 X 10-* 


N 




9.65 X 10-5 


9.65 X 10-5 


O 




8.55 X 10-* 


8.55 X 10-* 


Ne 




3.84 X 10-5 


3.84 X 10-5 


Mg 




4.17 X 10-5 


4.17 X 10-5 


Si 




4.94 X 10-5 


4.94 X 10-5 


S 




1.64 X 10-5 


1.64 X 10-5 


Fe 


3.7 X 10-5 


3.29 X 10-5 


6.58 X 10-5 


X 


0.695 


0.712 


0.711 


Y 


0.304 


0.270 


0.269 


Z 


2.066 X 10-2 


2.624 X 10-2 


2.808 X 10-2 



sidered the model of solar composition with doubled number 
abundance of iron named as Comp. III. 

In this paper, we have calculated our canonical model of 
the disk with chemical composition Comp. I, and the power- 
law irradiation with spectral index ax = 0.9, and X-ray 
luminosity Lx = 10'*'' erg s-^. For chemical composition, 
Comp. I, we compared this canonical model to the case of 
power-law with the same spectral index, but different lumi- 
nosities: Lx ~ 10*2 erg and Lx = 10** erg , and to 
the case of power-law with the same luminosity but different 
spectral indices: ax = 0.6 and ax ~ 1.2. We also compared 
the structure and spectrum of our canonical model to models 
with the luminosity Lx ~ W*^ erg s^^ , and spectral index 
ax = 0.9, but with other chemical compositions: Comp. II 
and Comp. III. In all computed models the power-law radi- 
ation extends between hvmin = 0.1 and hvmax = 100 keV. 



3 EQUATIONS OF AN IRRADIATED 
ATMOSPHERE 

The set of equations presented and solved here represents a 
direct extension of the equations presented in MR04. Now 
we added terms which correspond to the ffuorescent Ka and 
Kf} lines of iron, adjusted to the accretion disk geometry and 
assumed that the external irradiation from a point X-ray 
lamp has the power-law spectral distribution. 
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3.1 The equation of transfer 



3.1.3 Compton scattering emission 



The equation of transfer for tiie specific intensity I^, at fre- 
quency V was written in plane-parallel geometry on the ge- 
ometrical depth scale z 



(1) 



M = Ji/ - (liiy + )U , 

paz 

where j^,, and denote frequency dependent emission, 
absorption and scattering coefficients for 1 gram, respec- 
tively. In this paper we use the LTE (local thermodynamic 
equilibrium) absorption k^, whereas coefficients of emission 
jv and scattering a,^ include nouLTE terms. 

Eq. [1] can be written on the monochromatic optical 
depth, with dr^ = —{ku + (jy)pdz, then 

^^=^--^^=^^-^- (2) 



where 5*^ is the frequency dependent source function. 

Emission coefficient jv is the sum of three terms, jv ~ 
3v^ + jf^ + jv\ which represent thermal emission, Compton 
scattering emission and the emission in iron fluorescent lines, 
respectively. 



3.1.1 Thermal emission 

Coefficient of thermal emission in LTE is proportional to 
the Planck function 

fj" = . (3) 

The coefficient of true absorption the sum of bound- 

free absorption from numerous levels of atoms and ions for 
all elements, plus free-free absorption from all ions. 

We also included absorption of 4 lowest lines of funda- 
mental series of helium-like iron and of similar 4 lowest lines 
of hydrogen like iron, all formed in LTE by assumption. 

3.1.2 External irradiation 

We defined two angle-averaged variables, which were derived 
from the specific intensity 7J^* of the external irradiation 
penetrating the atmosphere to the depth 



(4) 



Vl,{t„) = I 7*^*(cj,r^) exp (— r^Z/i^) dcj , (5) 
4^ In 



where p.^: stands for the cosine of zenithal angle of the run- 
ning direction ui. The external intensity from the point-like 
lamp is emitted in the form of power-law with spectral index 
ax: 



Text A —Oiv f ^ \ { ^min \ 

1,^ = Av ^ exp exp , 

V Umax J \ V J 

and normalized to the luminosity of the source, Lx '■ 

Lx 



A 



4nrf 



(6) 



(7) 



The distance from an X-ray source depends on the ring ra- 
dius r in the casual relation rf = + r^. The luminosity, 
spectral index, and cut-off limits of irradiating spectrum are 
free parameters of our model, all of them described in Sec. (2] 



Expressions for the Compton scattering terms and emission 
coefficient in an irradiated stellar atmosphere were derived 
in MR04 and are valid in a disk atmosphere without any 
changes 



(8) 



In the above equation variable Uv denotes the angle- 
averaged intensity of the external irradiation (Eq. [Jjl . 

Compton scatte ring cross sectio ns were computed fol- 
lowing the paper bv iGuilbertI (|l98ll ). Functions $i and $2 
are properly weighted angle- averaged Compton redistribu- 
tion functions for photons both incoming or outgoing of fre- 
quency V after scattering in thermal electron gas, cf. MR04. 



3.1.4 Iron fluorescence emission 

Fluorescence of low-ionized iron gas was approximated by 
two emission lines, Ka and Kfj. Therefore, 



(9) 



where _E£' and E^^ denote the integrated intensity of Ka 
and Kp emission lines, respectively. Frequency dependent 
variables and tp^ define profiles of fluorescent lines, both 
normalized to unity. 

Energy emitted in and Kjj lines was derived by ab- 
sorption of hard continuum X-rays from the radiation field 
penetrating the disk atmosphere. Hard X-ray photon inter- 
acting with neutral or low-ionized iron most probably ionize 
and remove electron from the innermost K shell, and then 
remaining hole is filled by another electron falling from L 
shell [Ka) or M shell (Kp transition). The integrated emis- 
sivity for Fe fluorescent line is given by: 



hv 



■{Jv + Uv)dv, 



(10) 



where we set Y = 0.34 (jKrause fc Olived 1 19791 1. Variable 
a'J°" is the bound-free absorption coefficient for ionization 
from K shell of iron counted for 1 atom. Similar expression 
holds for the integrated intensity E^^ . 

The fluorescence yield, Y , defines the fraction of the 
energy of hard X-rays absorbed by iron atoms which was 
reemitted as photons in Ka line. The remaining energy of 
absorbed X-rays, 1 — F, was spent for the release of Auger 
electrons. 

The iron Ka fluorescent line is a doublet line, and such 
a structure was reproduced by our code. We set the central 
energies for Ka^ and Ka^ lines to 6.404 keV and 6.391 keV, 
respectively. Natural widths (FWHM) of both lines, 2.7 eV 
Kg^ ) and 3.3 eV (Ka^i) were taken from iKrause fc Oliverl 
I979I ). The Kp line was approximated by a singlet line at 



the energy 7.057 keV. We set the natural width of the line 
to 2.5 eV arbitrarily. 

Integrated emissivity Ea was divided in the proportion 
2:1 bet ween Kg^ an d Ka.2 components of the doublet. Fol- 
lowing [Bask^ l| 19781 ) we quite arbitrarily assumed that the 
integrated intensity _E{' = 0.1 Ea. Opacity proflles of all 
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three lines were set to Voigt functions with depth-dependent 
parameters describing natural and thermal broadening. 



3.2 Equation of state 

We used here the LTE equation of state, where the number 
of free electrons and all ionization populations of iron was 
determined by the Saha ionization equations. 

In the case of an atmosphere irradiated by strong non- 
local radiation field the assumption of LTE yields inaccurate 
ionization and excitation populations. For instance we get 
inaccurate populations of iron atoms with various configu- 
rations of electrons in the outermost M and N shells. 

This is the reason that we did not discuss in our paper 
intensities for the spectral features of energies below the 
soft X-ray excess in the outgoing spectra, where the nonLTE 
configuration of outermost bound electrons in various atoms 
is of vital importance. 

However, we believe that the population of iron atoms 
with filled innermost K and L shells is always equal to the 
total number of iron atoms with 8 or more bound electrons, 
and does not depend on the LTE/nonLTE details of elec- 
tron populations in M and N shells. As the result we believe 
that have obtained reasonable intensities for both fluores- 
cent lines Ka and Kp of low-ionized iron. The lines were 
considered strictly in nonLTE. 



3.3 Equations of hydrostatic and radiative 
equilibrium 

3.3.1 Hydrostatic equilibrium 

We solve the static equation of equilibrium for gas pressure 

dPgas Q dPrad 



dr 



(K + a) 



dr 



(11) 



which takes into account gradient of radiation pressure Pr 
at the level r. Expanded form of the equation is 



HP 

gas 

dr 



(k + a)std 



— / rj^iHi, - K) du , 



(12) 



where 77^ = (k^ -|- o"i/)/(k -I- a)std. Eq. 1121 neglects possible 
vertical motion of matter which can be caused by the differ- 
ential rotation of the disk. 

In addition we assumed that the vertical gravitational 
acceleration in the atmosphere of a given ring is constant i.e. 
its dependence on the height above the disk equatorial plane 
was neglected. The value of acceleration was computed from 
the solution of the disk model presented in Sec. [21 



3.3.2 Radiative equilibrium 

Constraint of radiative equilibrium requires that 



Hv{ti,) dv 



Vv(t,j) dv 



47r 



(13) 



at each depth level of in the atmosphere. Radiative con- 
stant aR — 5.66961 x 10^^ (cgs units). In this paper, we 
assumed that whole flux from an accretion disk was gener- 
ated below an atmosphere, and we do not solve local energy 
generation via viscosity. 



The above equation subsequently was used to determine 
distribution of temperature in the disk atmosphere at a given 
ring. Transformation of Eq.[T3]to the expression for tempera- 
ture corrections ST was done by the linearization technique. 
The equation of transfer including temperature corrections 
was solved using t he technique of variable Eddington fac- 
tors (|Mihalaslll978l ). Details of the full set of equations used 
in our computations and the method of solution by partial 
linearization are presented in MR04. 



4 COMPUTATIONAL RESULTS 

Our computer code allowed us to compute the structure of 
disk atmospheres over very large range of electron scattering 
optical depth starting from Tes = 10~* up to res = 10^. Fur- 
thermore, we were able to reproduce the overall continuum 
spectrum from deep infrared of 0.4 eV up to hard X-rays 
of 400 keV. We present our spectra as an energy dependent 
outgoing specific intensities, 1^, which are suitable for disk 
geometry. We reject presentation of monochromatic fluxes 
since they are relevant only to geometry of a spherical star. 



4.1 Temperature structure 

Fig. [1] represents the temperature structure of various ir- 
radiated models. The canonical model for several rings is 
presented in panel (a). In panel (b) our canonical model is 
compared with models of the same spectral index, but with 
different X-ray luminosities, while in panel (c) with models 
of the same X-ray luminosity but with different spectral in- 
dices. Panel (d) shows the structure of the disk irradiated 
by the power-law with ax = 0.9 and luminosity Lx ~ lO^'' 
erg s"'^ but for various chemical compositions. In panels (b), 
(c), (d) we show temperature structure only for the inner- 
most ring number 1. 

In all computed cases we note the existence of very 
hot outermost at mospheric layer which is caused by ex- 
ternal irradiation ([N ayakshin ct al.' '2000'; Balla ntvne et all 
,20011 : iHubenv et aU .2001: Roz anska ct al, 2002). The effec- 
tive temperature of all rings resulted from the accretion rate 
and non-zero viscosity, is slightly less or equal to 10^ K. How- 
ever, the temperature of the outermost hot skin approaches 
3 X lO'^ K. For deeper layers external irradiation practically 
does not change the temperature of the gas. 

The surface temperature is determined by the balance 
between Compton heating caused by hard X-ray photons 
from the external illumination and inverse Compton cooling 
by the low-energy reprocessed radiation and thermal radia- 
tion from the disk below. Hence the actual surface temper- 
ature is lower than the Compton temperature due to the 
illumination alone. 

In our canonical model, the temperature of the hot skin 
is almost the same for all rings (panel (a)) and apparently 
it does not change with the distance from the central black 
hole for a given power-law hardness and the luminosity of 
the external irradiation. We demonstrate, that the surface 
temperature depends on the X-ray luminosity (panel (b)) 
and the hardness of spectrum (panel (c)). For softer power- 
law, ax = 1.2, Compton heated skin has lower temperature, 
but transition from hot to colder layer occurs deeper in the 
atmosphere (long dashed line in panel (c)). 
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Figure 2. Ionization structure of three main elements for the 
model with solar chemical composition (Comp.II see Tab. [2]l. 
Panel (a) shows ionization fraction of carbon ions, panel (b) - 
oxygen ions, and panel (c) - iron ions. 



Figure 1. Comparison of temperature structure for different 
models. In panel (a) our canonical model is presented for ring 1 - 
solid line, ring 4 - short-dashed line, ring 5 - dotted line, and ring 
8 - long-dashed line. In panel (b) canonical model (solid line), 
is compared with models of different luminosities: Lx = 10^'* 
erg s~^ - long-dashed line, Lx = 10"^^ erg s~^ - short dashed 
line, and dotted line - model with no irradiation. The same canon- 
ical model (solid line) in panel (c) is compared with the model 
of ax = 0.6 - short dashed line, and ax = 1-2 - long dashed 
line. Panel (d) shows the canonical model (solid line), compared 
with the structures of the atmosphere illuminated by the same 
radiation: Lx = 10*'^ erg s~^, ax = 0.9, but for Comp.II - short 
dashed line, and Comp. Ill - long dashed line. Only first rings are 
presented for all models in panels (b), (c) and (d). 



Table 3. Description of lines in Figs. [3} [8] 



type of a line 


cos('t) 


i 


solid black 


0.9801 


11.4° 


short-long dashed 


0.8983 


26.1° 


long dashed dotted 


0.7628 


40.3° 


short dashed dotted 


0.5917 


53.7° 


dotted 


0.4083 


65.9° 


short dashed 


0.2372 


76.3° 


long dashed 


0.1017 


84.2° 


solid red 


0.0199 


88.9° 



The transition between hot and cold zone is very steep, 
geometrically thin, and depends on the initial conditions of 
a particular ring. For the first ring the transition occurs 
deeper in an atmosphere than for the outermost eight ring 
(panel (a)). The shape of the transition is caused by ion- 
ization on heavy elements and sometimes possesses small 
horizontal branch, which is consistent with results from pre- 
vious photoionization computations (Navaksh in fc Kallm an 
I2OO j ). The location of such horizontal branch changes for 
diflFerent chemical compositions as shown in panel (d). 

The temperature profile of an irradiated disk causes the 
ionization structure of the atmosphere to exhibit two dis- 



tinctly separated regions. The surface high temperature zone 
is just the region where heavy elements are in highly ionized 
states. Below that all elements change into low ionization 
states. We have illustrated this behaviour in Fig. |2l pre- 
pared for the ring 1, irradiated by power-law with ax = 0.9 
and luminosity Lx = 10*^ erg s~^, with solar chemical com- 
position, Comp. II, (see: Tab. [2]). The ion number denotes 
the number of missing electrons in particular atom. In the 
outermost hot layer, oxygen and carbon are completely ion- 
ized, while the most abundant iron, Fe24, is in the form of 
helium-like. 
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Non — irradiated 
disk 




-2 2 

log Energy (keV) 

Figure 3. Outgoing spectra from the disk integrated over radii for 
eight different aspect angles (see Tab.|3]l. Red soUd line represents 
edge on disk, while black solid line represents disk seen face on. 
Panel (a) shows spectra of the disk with no irradiation. Panels 
(b), (c), (d) show models with the same ax = 0.9, but for three 
different luminosities 10*'^, 10*^ and 10^'* erg s~^, respectively. 
The most important ions producing the key spectral features have 
been indicated. 



4.2 Outgoing intensity spectra 

Figures |3]-|8]present outgoing intensity spectra of the whole 
accretion disk at various aspect angles. The spectra were 
obtained by integration over eight rings located at diflerent 
radii, see Tab. [1] Particular lines in Figs. |3] - |8] correspond 
to different angles between the direction to the observer and 
the normal to the disk. Exact values of those eight angles 
and their cosines are given in Table [S] In further discus- 
sion we draw attention of the reader to the extreme angles: 
solid black line represents almost vertical direction (face-on 
aspect), whereas the solid red line represents almost hori- 
zontal direction (edge-on aspect). 

The spectrum of our canonical model is compared with 
the model of the same spectral index but different X-ray 
luminosities. Figs. |3l and|6] and of the same X-ray luminos- 
ity but different spectral indices Fig. 3] and [T] Spectra for 
various chemical compositions and the same irradiation are 
plotted in Figs. [5] and [8] 
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Figure 4. Outgoing spectra from the disk integrated over radii 
for eight different aspect angles (see Tab.fSjl. Red solid line repre- 
sents edge on disk, while black solid line represents disk seen face 
on. Panels (a), (b), (c) show models with the same Lx = Iff*^ 
erg , but for three different spectral indices 1.2, 0.9 and 0.6, 
respectively. 



Panel (a) in Fig. |3] presents intensity spectra of the disk 
with no irradiation. Comparison to panels (b)-(d) shows, 
that all emerging radiation above above 0.2 keV originates 
from the disk atmosphere as the result of an external hard 
X-ray illumination. 

In all irradiated models outgoing spectra show the ef- 
fect, that the intensity of radiation emerging almost hori- 
zontally to the disk {i = 88.9°) is greater than the intensity 
emerging almost vertically {i — 11.4°). This is valid in the 
X-ray energy band from about 0.1 keV up to 100 keV, and in 
the infrared band. Such an effect of limb-brightening should 
be attributed to the inversion of temperature in the irradi- 
ated disk atmospheres, see Fig. [T] 

All broad-band spectra show the deficit of hard X-rays, 
which were absorbed in the disk and suffered continuum 
Compton downscattering. This caused that the energy of 
hard X-rays was redistributed to lower energies. Moreover, 
there exists a spectral bump below 1 keV, which is similar 
to the well known soft X-ray excess. Such bump is more 
prominent for the disk seen edge on (red solid line) in all 
models. 



5 IRON LINE COMPLEX 

External irradiation generates two iron fluorescent lines: Ka 
doublet at 6.391 and 6.404 keV, and singlet at 7.057 keV. 
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5.1 Chemical composition of the accretion disk 

Presence of iron atoms in the gas of Comp. f ensures gen- 
eration of iron Ka and iron fluorescent lines and their 
Compton Shoulders (Tab.|3]flrst five rows). Both fluorescent 
lines originate in deeper (colder) atmospheric layers, where 
iron atoms are only weakly ionized. However, if the actual 
chemical composition is rich of various heavy elements, like 
C, O, Si, S etc., then the appearance and equivalent widths 
of both fluorescent lines can be reduced due to K-edge ab- 
sorption of X-ray photons by all those elements. 

To check this we have computed disk models for our 
canonical X-ray luminosity Lx ~ 10*'' erg s~^ and spec- 
tral index ax ~ 0.9 assuming other chemical compositions 
(Comp. II and Comp. Ill) with the excess of helium and the 
solar chemical abundance of eight heavier elements: C, N, 
O, Ne, Mg, Si, S and Fe (see Table. 0. 

Table. |4] shows that the equivalent widths of the Comp- 
ton Shoulder has the highest value in models of Comp. I, 
both for Ka and Kp lines. Addition of heavy elements other 
than iron (Comp. II, and Comp. Ill) and their K-edge ab- 
sorption, significantly decreases EW of Compton Shoulder 
(compare rows 2, 6 and 7 of Table. |4]). The same holds for Ka 
and Kp line cores seen edge-on. But for the disk seen face-on 
{i = 11.4°), the situation for line cores is opposite. EWs of 
both line cores increase with addition of heavy elements. 



Figure 5. Outgoing spectra from the disk integrated over radii for 
eight different aspect angles (see Tab.|3]l. Red solid line represents 
edge on disk, while black solid line represents disk seen face on. 
Panels (a), (b), (c) show models with the same ax = 0.9, and 
Lx = 10'*'' erg s~^, but for three different chemical compositions: 
Comp. I, Comp. II, Comp. Ill respectively, given in Tab. [2] 



Furthermore, due to Compton scattering of the line pho- 
tons in the irradiated atmosphere photons from line core 
are shifted toward lower energies creating red wing, which is 
called Compton shoulder. All the above features are simul- 
taneously reproduced in our model calculations. 

Figs, ini [3 [H] present integrated intensity spectra just 
around the fluorescent Ka doublet (left panels), Compton 
shoulder (middle panels), and fluorescent Kfs singlet (right 
panels) for all computed spectra. Unfortunately, Compton 
shoulder of Kf3 line is very weak in our models. 

Limb-brightening is clearly seen in the region of iron 
line complex. The Compton shoulder is most prominent in 
the vertical direction (black solid line), while it almost dis- 
appears for the disk seen edge on (red solid line). 

The amount of energy emitted by iron line features de- 
pends on the shape and luminosity of irradiating contin- 
uum and on its spectral index. We have computed equivalent 
widths of modeled iron line components, and collected them 
in Table m for two extreme aspect angles: face on (i = 11.4°) 
and edge on, i — 88.9°. Features are stronger for models 
with lower luminosity and harder spectral index. For such 
cases, also Ka and Kp lines seen edge on are more promi- 
nent. The energy emitted in Compton shoulder is about ten 
and more times lower than the energy emitted in Ka line. 



6 DISCUSSION 

Calculations of an accretion disk atmosphere irradiated by 
external hard X-rays are very time consuming and some- 
times converge with difficulty. In our paper we have com- 
puted the disk up to the distance from the central black 
hole equal ISRschw ■ Model atmospheres of farther rings did 
not converge with our computer code. Going farther from 
the black hole, the gravity in the disk atmosphere becomes 
lower, and irradiation disturbs hydrostatic equilibrium due 
to strong radiative pressure. This problem should be solved 
using dynamical calculations, but none of recently avail- 
able hydrodynamic codes contains such complicated radia- 
tive transfer as presented in this paper. 

Our treatment of Compton scattering allowed us to 
compute the structure of the Compton heated skin in an 
accretion disk atmosphere toge ther with the transit ion layer 
down to the deep cold zone l|Hubenvet al.ll200ll ). Struc- 
ture of the transition layer depends on heavy elements abun- 
dance. We have shown that the position of a small horizontal 
branch around 3 x 10^ K changes with chemical composition 
of an irradiated gas. 

In the soft X-ray domain of our model spectra we ob- 
serve huge increase of speciflc intensity due to Compton re- 
distribution of hard X-ray photons toward lower energies. 
We claim here, that those bumps give impact to soft X-ray 
excess frequently observed in many objects. Such excess in 
our models is higher for higher luminosity of X-ray source 
and for softer (i.e. higher) power-law index. 

All spectra discussed here were computed in the frame 
of the accretion disk, so we did not consider line broaden- 
ing and distortion caused by Keplerian motion of gas around 
black hole. Our goal was to demonstrate the existence of the 
red wing due to Compton scattering of fluorescent line pho- 
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Table 4. Equivalent widths of modeled iron line components: Ka , its Compton Shoulder (CSa), and Kp, and its Compton Shoulder 
(CSp). All values are given in eV for two extreme aspect angles: i = 11.4° - face-on, and i = 88.9° - edge-on. 
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Figure 6. Iron line complex in the disk spectra integrated over radii for 8 different aspect angles (see Table (Sj. Panels (a), (b),and (c) 
show model spectra with the same ax = 0.9, but for three different luminosities lO^^, 10*3 and 10** erg s ^ , respectively. Iron Ka 
doublet is presented in the left panels, whereas Compton shoulder of this line is shown in the middle panels. Right panels show Kfj line 
with Compton shoulders of a marginal height. 



tons on a colder electron gas. Narrow Ka lines recently ob- 
served by XMM Chandra and Suzaku sat ellites have equiva - 
lent width 75 eV in ca se of MCG-H8-11-11 l|Matt et al.ll2006h . 
165 eV in NGC 2992 (lYaaoob et al.ll2007l) and 190 eV in the 
famous object MCG-6-30-15 l|Miniutti et all [2007V Those 
values of EWs are reproduced by our model calculations, 
see Table |4] The resolution of actually working satellites is 
unable to resolve iron Ka double profile. 

Similar comparison with Kp line s hows excellent agree- 
ment of our models with observations. lYagoob et all (|2007l ) 
reported detection of Kg iron lin e in N GC 2992 with an EW 
of about 27 eV, and iMatt et all (|2006l ') estimated an upper 
limit of 10 eV in MCG-h8-ll-ll. 

There are not many detections of Compton shoulder in 
AGN. Most probably, because it is very difficult to distin- 
guish between red wing produced by Compton scattering 
and that caused by relativistic disk motions. Ne v erthel ess, 
EW of Compton shoulder fitted by iMatt et al.l (120061 ) in 
MCG+8-11-11 equals 17 ± 10 eV. Such a value agrees with 
values obtain in our computations for the disk seen face on, 



which contains only hydrogen, helium and iron (Comp. I). 
EW and the profile of Compton shoulder presented in this 
paper were computed on the bases o f deta iled model atmo- 
sphere calculations, whereas iMattI l|2002h estimated them 
from models of irradiated "wall" of constant density. 

Our models of broad-band intensity spectra from the ir- 
radiated accretion disk are interesting also in other spectral 
domains, for example in infrared and optical bands. They 
can be immediately used to study other important spectral 
features like hydrogen Lyman edge which looks different de- 
pending on aspect angle. We address this issue for future 
consideration. 



7 SUMMARY 

In this paper we presented intensity spectra emerging from 
a sample accretion disk around a supermassive black hole in 
AGN. The disk was irradiated by hard X-rays of the power- 
law spectral shape originating from a point-like source lo- 
cated above the inner disk. We compared results for the 
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Figure 7. Iron line complex in the disk spectra integrated over radii for 8 different aspect angles (see Table IS}. Panels (a), (b),and (c) 
show model spectra with the same Lx = W^^ erg s~^, but for three different spectral indices 1.2, 0.9 and 0.6, respectively. Iron Ka 
doublet is presented in the left panels, whereas Compton shoulder of this line is shown in the middle panels. Right panels show K/j line 
with very low red-wing Compton shoulders. 
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Figure 8. Iron line complex in the disk spectra integrated over radii for 8 different aspect angles (see Table [Sjl- Panels (a), (b),and (c) 
show model spectra with the same Lx = 10*'^ erg s~^, and spectral index ax = 0.9, but for three different chemical compositions, 
Comp. I, Comp. II, and Comp. Ill respectively. Iron Ka doublet is presented in the left panels, whereas Compton shoulder of this line 
is shown in the middle panels. Right panels show enlarged Kjj lines with their weak Compton shoulders. 



cases when irradiation luminosity, power-law spectral index, 
and chemical composition change. 

Reprocessing of external hard X-rays causes generation 
of a very hot external layer in the disk atmosphere where 
temperature rises to millions K, whereas the structure of 
deeper atmospheric layers is not affected at all. 

We numerically reproduced redistribution of external 



X-rays both by thermal emission from the heated disk skin 
and by Compton downscattering in deeper and colder lay- 
ers of the disk atmosphere. As the result outgoing intensity 
spectra of the disk exhibit the following important observ- 
able properties. 

First, external irradiation generates two fluorescent 
hues of iron, Ka and K/s. Compton scattering of the line 
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radiation on colder electron gas (below the hot disk skin) 
generates the Compton shoulder on the red side of both Kc, 
and Ki3 fluorescent lines, nevertheless shoulder is very weak 
for K/s line. Compton shoulders m odeled here differs from 
the Compton shoulder presented bv lMatd l|2002l ). when the 
latter resulted from X-ray reflection off the constant den- 
sity slab. Equivalent widths of iron lin e complex agree with 
recent ob servations (|Matt et all |2006| : lYagoob et al.l l2007l: 
iMiniutti et al. 200?! '). 

Second, Compton scattering of the external X-rays 
causes build up of the soft X-ray excess below 1 keV. This 
excess is higher for the disk observed face on, but still too 
low to explain observations. 

Third, our results clearly demonstrate the effect of limb 
brightening of the disk radiation mostly in the X-ray do- 
main. This result can put new constrains on the geometry 
of the X-ray source and the reflecting disk atmosphere. We 
address this issue for further consideration. 
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